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rH . Abstract 

The coupled Einstein-Dirac-Maxwell equations are considered for a static, 
Q>^ ■ spherically symmetric system of two fermions in a singlet spinor state. Stable 

soliton-like solutions are shown to exist, and we discuss the regularizing effect 
of gravity from a Feynman diagram point of view. 
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cn . There are some interesting effects that result when one couples gravity, as ex- 

pressed through Einstein's equations, to other fundamental force fields. The most 
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O I striking example of this is the discovery by Bartnik and McKinnon (BM) [|T| of 

soliton-like solutions for the SU{2) Einstein- Yang/Mills (EYM) equations. This 
came as a surprise because several results for related systems lead to the conjec- 
Q^. ture that such solutions cannot exist . For example, neither the vacuum Einstein 

5_i ! equations, nor the Yang/Mills equations admit non-trivial static, globally regular 

solutions. The existence of the BM solutions, however, is due to the coupling of 
different force fields, whereby the attractive gravitational force is balanced by the 

^ ■ repulsive Yang/Mills force. But this balance is quite delicate, so much so, that the 

EYM solutions are unstable to small perturbations [^]. 

We consider here a different type of coupling; namely gravity coupled both to 
quantum mechanical particles with spin (Dirac particles) and to an electromag- 
netic field. We analyze the resulting Einstein-Dirac-Maxwell (EDM) equations for 
a static, spherically symmetric system of two fermions in a singlet spinor state. We 
construct stable soliton-like solutions, and we discuss some of their properties for dif- 
ferent values of the electromagnetic coupling constant. Our aim is to show that the 
gravitational interaction has a certain "regularizing" effect; this will be explained 
from a Feynman diagram point of view. 

Before going further, we point out that throughout this essay, we are considering 
only classical fields coupled to the Dirac equation; in particular, we consider neither 
quantized bosonic fields, nor second quantization of the Dirac field. 

*This essay received an "honorable mention" in the 1999 Essay Competition of the Gravity 
Research Foundation - Ed. 
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^Research supported in part by the NSF, Grant No. 33-585-7510-2-30. 



The EDM equations for a system of n Dirac particles have the following general 
form, 
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where j^ = eX]a=i ^iG^^^a is the Dirac current, and T- is the sum of the energy 
momentum tensor of the Dirac particles and of the Maxwell stress-energy tensor. 
Here Kj is the Ricci tensor, Fjk is the electromagnetic (em) field tensor, and "^a are 
the 4-component Dirac wave functions describing particles of mass m and charge 
e. The Dirac matrices G^ are (4 x 4) matrices, which are related to the Lorentzian 
metric via the anti-commutation relations ^{G^{x), G^{x)} = g^'^lx)!. Dj is the spin 
connection; it is composed of the ChristofFel symbols of the metric connection and 
of the electromagnetic potential (see [H). The EDM equations only make physical 
sense if the wave functions \l'a are properly normalized. For this, we choose a space- 
like hypersurface 7i with a future-directed normal vector field z/ and demand that 



/ '^aG^'^b^jd^i = 6ab , a, 6=1 

Jn 



where dfi is the invariant measure on 7i, induced by the Lorentzian metric. 

The Dirac and Maxwell equations are linear in the wave functions and in the 
electromagnetic field, respectively. But the Dirac current and the energy-momentum 
tensor are quadratic in the "^a, the stress-energy tensor is quadratic in the electro- 
magnetic field, and the Einstein equations are nonlinear in the metric. Thus the 
EDM equations form a highly nonlinear system of differential equations, which in 
this generality is far too complicated to analyze. Therefore we consider a simple 
ansatz where the EDM equations reduce to a system of ordinary differential equa- 
tions; namely a static, spherically symmetric system of two Dirac particles in a 
singlet spinor state. For the wave functions, "static" means that their time de- 
pendence is a plane wave e~*'^*. The parameter tu > is the energy of the Dirac 
particles. Using a special rescaling technique [^ ^ , we construct numerical solutions 
of our EDM system. 

We now briefly mention some general properties of the solutions. All of our solu- 
tions are everywhere regular, asymptotically Minkowskian, and have finite (ADM) 
mass p. The Dirac wave functions decay exponentially at infinity, so that the 
Dirac particles are with a high probability confined to a neighborhood of the origin; 
this is due to the gravitational attraction. For large r, our solutions look similar 
to a Reissner-Nordstrom solution with charge 2e and mass p. Our solutions are 
parametrized by two physical parameters: the rest mass m of the Dirac particles 
and the relative strength of the electromagnetic and gravitational coupling (e/m)^. 
For a fixed value of these parameters, we find an infinite, discrete number of so- 
lutions which can be regarded as the ground state, the first excited state, etc. In 
Figure |I|, the mass-energy spectrum for the ground state; i.e., the binding energy 
m — u versus m, is plotted for different values of (e/m)^; these curves are spirals. 
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Figure 1: Binding Energy m — u oi the Fermions for (e/m)^ = (A), 0.7162 (B), 
0.9748 (C), 1 (D), and 1.0313 (E). 



The binding energy decreases for fixed m as {e/uif' increases; this means that the 
em repulsion weakens the binding. For {e/m) > 1, the electromagnetic repulsion is 
so strong that there are no solutions in the Newtonian limit. However, due to the 
nonlinearity of Einstein's equations, the EDM equations still admit solutions (spiral 
E). By combining linear perturbation theory with more subtle topological methods 
(Conley index, cf. 0), the stability of the solutions is analyzed |^. For (e/m)^ < 1, 
we find that all solutions on the "lower branch" of the spiral (i.e. on the curve from 
m = up to the maximal value of m) are stable, whereas all solutions on the "upper 
branch" are unstable. 

Our solutions of the EDM equations have some interesting features. In this essay, 
we shall focus on the discussion of what the form of the spiral curves in Figure |l] tells 
us about the nature of the gravitational interaction (see p, ^ for the discussion of 
other effects). If the EDM equations (|I]) are treated in a perturbation expansion (in 
the spirit of [§), we obtain a description of the interaction with Feynman diagrams. 
More precisely, we get all the Feynman diagrams of quantum field theory, with the 
only exception of the fermionic loop diagrams; these are absent because we do not 
consider pair creation/annihilation (the number of Dirac particles is fixed in our 
ansatz). Thus we can say that, although we consider only classical fields, many 
aspects of perturbative quantum field theory are already incorporated in our EDM 
equations. The first question arising from this point of view is why the divergences 
of the Feynman diagrams do not show up in our regular, well-behaved solutions. 
To analyze this problem in a simple example, we consider the loop diagram of the 
virtual emission and absorption of a photon in Figure |^. If evaluated in Minkowski 
space, this diagram involves the divergent integral over the product of two Green's 




Figure 2: A simple loop diagram 
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In our ansatz, we only consider a static electromagnetic potential, which means that 
the momentum q of the electromagnetic Green's function must be purely spatial, 
q= (0, g). Thus the process is described, instead of (0), by the three-dimensional 
integral 
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This integral is finite, as one sees by considering the behavior of the integral for large 
|g|^, taking into account spherical symmetry. This illustrates that the reason why 
our interaction is finite is because of our restriction to static, spherically symmetric 
solutions; note that it is unrelated to the fact that our fields are merely classical. 

We now come to our main point: the convergence of the perturbation expansion. 
A perturbation expansion in the gravitational coupling constant corresponds to a 
Taylor expansion of our system in the mass parameter m. We consider as an exam- 
ple a spiral of Figure |I]for small (fixed) (e/m)^. The just explained finiteness of all 
Feynman diagrams becomes manifest by the fact that the spiral is smooth near the 
origin. However, since it is not even one-to-one, the spiral is not analytic in m. This 
implies that the perturbation expansion cannot converge absolutely for large values 
of the coupling constant. This explicit non-convergence result is a peculiar effect of 
the gravitational interaction; e.g. the mass-energy spectrum for the Dirac-Maxwell 
system is analytic. Nevertheless, since the spiral curves are bounded, we see that 
the total interaction is finite even for strong coupling. This implies that there must 
he cancellations between Feynman diagrams of different order, such that the pertur- 
bation series is finite, although this series does not converge absolutely. Thus the 
gravitational interaction of the Dirac particles is more regular than what one would 
expect from a naive perturbation expansion. Actually, it has long been conjectured 
that the ultraviolet divergences of quantum field theory might disappear when grav- 
itation is included in a non-perturbative way (see [^ for a recent formulation of 
this conjecture). Our EDM system allows us to verify this conjecture explicitly in 
a simple example. Our results also show that effects like conditional summability 



of the perturbation series actually occur for the fundamental physical interactions 
with gravity. 
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